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ABSTRACT 

Recent STIS and FUSE observations have detected O VI absorption lines at low redshift that are not 
clearly associated with any galactic system. In this paper, we argue that these lines are due to metal 
enriched hot gas in the intergalactic medium. Using numerical simulations of a cosmological-constant 
dominated cosmology, combined with reasonable assumptions about the metallicity distribution, we show 
that the number density and internal characteristics of these lines are correctly predicted. We find that 
the O VI is primarily produced by coUisional ionization from gas at a few times 10'^ K for lines with 
equivalent widths larger than 40 mA, while weaker lines can also be produced by photo-ionization. The 
absorption occurs in diffuse gas in filaments at moderate overdensity (5 ~ 5 — 100). 



Subject headings: 



intergalactic medium - 
— methods: numerical 
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1. INTRODUCTION 



The census of baryons in the local universe (Fukugita, 
Hogan & Peebles 1998) indicates that a substantial frac- 
tion of the baryonic density predicted by primordial nu- 
cleosynthesis remains to be detected. Numerical simula- 
tions (Cen & Ostriker 1999a; Dave et al. 2001) predict 
that a large fraction of this is in the form of a moderately 
warm/hot 10^ < T < 10^ K) component (Warm-hot 
intergalactic medium, or WHIM). The WHIM is hard to 
detect directly both because of the difficulty in UV obser- 
vations as well as due to its relatively low intrinsic density. 
A promising direction is to look for its signature in absorp- 
tion lines from background quasars (see Hellsten, Gnedin 
& Miralda-Escude 1998; Perna & Loeb 1998; Fang et al. 
2001). 

Recently, a number of groups have reported detecting 
the Li-like O VI resonance doublet in far-UV QSO spec- 
tra at z < 2 (Bergeron et al. 1994; Buries & Tytler 1996; 
Lopez et al. 1999; Churchill & Charlton 1999, rbhOl). Sav- 
age, Tripp & Lu (1998) also detected an O VI absorber at 
z = 0.225 without any C IV or N V absorption, indicat- 
ing that it was either coUisionally ionized with a temper- 
ature T ~ 5 X 10^ or photo-ionized with a low density. 
Tripp & Savage (2000) report an O VI absorption sys- 
tem a.t z — 0.14, possibly associated with a galaxy group 
and Tripp, Savage & Jenkins (2000) describe five O VI ab- 
sorbers in the range 0.15 < z < 0.27. Adopting reasonable 
assumptions, these authors argued that the total number 
of baryons associated with these absorbers was significant, 
comparable to the combined mass in stars, cool gas and 
cluster gas. 

In this paper, we explore the connection between these 
two ideas in more detail. In particular, we posit that the 
O VI absorption arises, in large part, in coUisionally ion- 
ized gas that resides in the sheets, filaments and low den- 
sity clumps that are naturally produced by gravitational 



collapse in hierarchical cosmology. The gas is shock-heated 
to a range of temperatures, while O VI is only significantly 
produced in a relatively narrow range around 3 x 10^ K. 
This is near the peak of the cooling curve (e.g. Sutherland 
& Dopita 1993) and so places a constraint on the density 
if the gas is not to cool in the Hubble timescale. For gas 
with 10% solar metallicity at 3 x 10^ K, the density must 
be less than nn < 10~^ cm""^. Assuming an O VI ioniza- 
tion fraction of 0.15, this requires a path length of 100 kpc 
to generate a column density of Nqvi = 4 x 10^'^ cm~^, 
typical of observed values. While this is a relatively large 
path length, the Hubble expansion if undisturbed would 
be only ~ 7 km s~^, which is well within the observed 
line widths. Therefore, we expect that thermal broaden- 
ing should dominate and b ^ 20 km s~^ for most lines. As 
we will show; however, colder, photoionized gas can also 
give rise to O VI absorption. For these lines, which tend 
to have lower column densities, the predicted line width 
is considerably smaller. Very recently, a preprint by Cen 
et al. (2001) has appeared which examines this issue using 
different simulation techniques; where there is overlap, our 
results agree reasonably well. 

This paper is organized as follows. In section 2 we de- 
scribe the hydrodynamic simulation and how we synthesis 
the simulated spectra. Section 3 discusses the results and 
Section 3 summarizes our main conclusions. 

2. METHOD 

We use numerical simulations to generate predictions 
for the observed O VI distribution. In particular, we use 
a grid-based Adaptive Mesh Refinement (AMR) method 
(Bryan 1996; Norman & Bryan 1999) that provides both 
good shock capturing as well as relatively high resolution 
in dense regions. We simulate a cube 20 h^^ Mpc on a 
side with gas and dark matter mass resolution of 6 x 10^ 
and 5 x 10* Af©, respectively. The smallest grid cell (i.e. 
the best resolution) is 9.8 kpc and occurs in the densest 
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regions. We assume a cosmological model with a mat- 
ter density flo = 0.3, baryon density fif, = 0.04, Hubble 
constant h — 0.67 (in units of 100 km s~^Mpc~^) and cos- 
mological constant f^A — 0.7. The initial density field is 
drawn from an adiabatic CDM power spectrum as approx- 
imated by the formula of Eisenstein & Hu (1998). 

We include dark matter and gravity but not radiative 
cooling or energy/metal injection from SN-driven galactic 
winds. By comparing a set of cosmological hydrodynamic 
simulations, Dave et al. (2001) concluded that radiative 
cooling play a minor role in the evolution of WHIM gas at 
z < 3 due to its low density. However, these processes can 
play an important role in metal enrichment of the inter- 
galactic medium (IGM). For instance, FUSE has detected 
O VI absorption in the galactic wind from the prototypical 
dwarf galaxies (see, e.g., Heckman et al. 2001). Since we 
do not know the distribution of oxygen a priori in this sim- 
ulation, we adopt a density dependent metallicity model 
as derived in the simulations of Cen & Ostriker (1999b). 
While this is reasonable, we plan to revisit this issue in a 
future paper using a simulation which self-consistently in- 
cludes both radiative cooling and metal/energy injection. 

To analyze the results, we generate artificial spectra as 
described in Fang, Bryan & Canizares (2001) and Zhang 
et al. (1997), both for the 1032 A O VI line and the 1216 
A HI Lya line. To compute the f{OVI/0) and f{HI/H) 
fractions, we adopt two models to investigate the effect 
of coUisional ionization and photo-ionization mechanisms. 
In model A, we use collisional ionization only, with ion- 
ization fractions from Mazzotta et al. (1998). In model 
B, we use Cloudy (version 90.04; Ferland et al. 1998) to 
generate a grid of temperatures (10^ < T < 10'' K) and 
densities (10~® < n(cm~^) < 10~^). We included a back- 
ground radiation field as computed by Haardt & Madau 
(1996) at z = due to the observed quasar distribution. 
We adopt a mean specific intensity at the Lyman limit of 
= 2 X 10"^'^ ergs s~"'^Hz~-'^sr~^. Since the simulation 
did not include radiative heating/cooling, the unshockcd 
IGM gas has a temperature which is unrealistically low. 
For these low temperature (and generally low density) re- 
gions, we adopt the relation between temperature and den- 
sity that has been found in simulations of the Lya forest: 
T = To{l + 5)^-1 where Tq - 5, 000 K, and 7 = 1.4 (see, 
e.g., Zhang et al. 1997; Ricotti, Gnedin, & ShuU 2000; 
Dave & Tripp 2001 for observational support). Here over- 
density 5 = pb/ {pb) — l and {pb) is the mean baryon density 
of the universe. To fit the simulated spectra we adopt the 
method developed by Zhang et al. (1997), which should 
work well for these relatively isolated lines. 

3. RESULTS 

Figure 1 displays the projected baryon overdensity 
(panel a), O VI column density (panel b, based on model 
B), temperature (10"* - 10^-^ K in panel c; 10^ - 10^ K in 
panel d). Comparing panel b with panel a, we find that 
most of the high column density O VI lines come from 
regions with overdensities oi 5 ^ 5 — 100. Collisional ion- 
ization is dominant in many of these regions, given their 
temperatures in the range 10^ — 10® K (panel d), but some 
systems are clearly photo-ionized, particularly at the low 
density end. Note that big clusters and groups are too hot 
to produce O VI. 
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Fig. 1. — Here we show various quantities from the simulation 
projected along a 5 ?t~^Mpc line-of-sight; the images are 20 
Mpc on a side. Panel a shows the m ean baryon overdensity, O VI 
column density (panel b, in units of cm~'^), temperature (10^ — 10^'^ 
K in panel c; 10^ — 10^ K in panel d). Most high column density 
O VI comes from regions with overdensities of 5 ~ 5 — 100 (panel 
a) and temperatures between 10^ — 10^ K (panel d). 

In Figure 2 we show the properties of the simulated re- 
gion along one random linc-of-sight (LOS). From top to 
bottom we display the baryon number density, peculiar 
velocity, temperature, O VI number density and transmis- 
sion spectrum {F = exp(— t), where r is the optical depth) 
for O VI in model A (solid line) , model B (dashed line) and 
for H I (dot-dashed line). Model B has been shifted by 5% 
for visual clarity. Note that at high densities and warm 
temperatures, both models agree well, while for low tem- 
perature regions, model B (which includes both photo- and 
collisional-ionization) , produces substantially more O VI. 

Figure 2e shows that both models generate a strong O VI 
absorption line at around z ~ 0.0019 along this LOS. No- 
tice that no strong peak appears in the n and novi dis- 
tributions at the corresponding redshift. This absorption 
line is actually produced by the novi peak at z ~ 0.0044. 
It is the peculiar velocity which shifts the line to a lower 
redshift: a peculiar velocity of ^ 750 km s^^ corresponds 
to a shift of Az ^ 0.0025. Parenthetically, we note that 
this means the total number of lines may be underesti- 
mated because the peculiar velocities can shift lines out 
of the simulated regions; however, we estimate that this 
would only change our predicted number density by about 
10%. Fitting the observed line in Figure 2e gives an equiv- 
alent width (EW) of ~ 185 mA and a Doppler 5-parameter 
of ~ 21 km s~^, which is roughly consistent with thermal 
broadening (6 = 0.129(r/A)i/2 km s~^ A is the atomic 
mass number). 
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Tripp & Savage (2000) suggested that a broad, hot H I 
Lya component may coexist with the O VI absorbers and 
such a broad component could be easy to miss with current 
observations. To demonstrate this, they showed that a fit 
to the miilticomponcnt Lya profile including a broad com- 
ponent (due to H I coexisting with the O VI in hot gas) is 
indistinguishable from a fit which does not contain such a 
broad component, at the signal-to-noise level of their data. 
In Figure 2e we display the H I Lya spectrum (the dot- 
dashed line) from our simulation. A broad component is 
clearly identified at the position of the strongest O VI ab- 
sorption line. Spectral fitting shows this line has a column 
density of 3.2 x lO^'' cm^^ and a Doppler ^-parameter of 
~ 76 km s^^. Both parameters are consistent with those 
from the missing H I Lya absorption line in Tripp & Sav- 
age (2000). The Doppler ?>-paramter is also consistent with 
thermal broadening by a hot gas of ~ 3 — 4 x 10^ K. ^ 
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Fig. 2. — Physical properties along one random LOS. (a) baryon 
number density; (b) pec uliar velocity; (c) temperature; (d) O VI 
number density (model A: solid line; model B: dashed line); (e) 
synthesized spectra for O VI and H I (dot-das hed line). In panel 
(d) and (e) model B is shifted down by 5% for visual clarity. The 
dashed line in panel (a) indicates the mean baryon density. The 
lines at ~ 0.0019 are produced by gas at 2 ~ 0.0044. 

By applying the spectral fit procedure to a total of 3, 000 
LOS we obtain the cumulative number per unit redshift 
[dN/dz) as a function of O VI equivalent width. Fig- 
ure 3 displays the dN/dz for both models. We also plot 
two observational results from Tripp & Savage (2000) for 
PG 0953-F415 and Tripp, Savage & Jenkins (2000) for 
E 1821+643. In general we find that the simulation and 
observations agree reasonably well. We notice that for 
EW > 40 mA, the O VI absorption line number per unit 
redshift from model B is less than twice the line number 
from model A. This implies that above EW = 40 mA 
the lines from coUisionally ionized gas begin to outnumber 
lines from photo-ionized gas. At EW > 80 mA, both mod- 



els give consistent distributions. This means that those 
strong lines come from high density and high tempera- 
ture regions, where collisional ionization dominates. For 
EW < 40 mA, lines from photo- ionized gas outnumber 
lines from coUisionally ionized gas, implying that lines pro- 
duced by photo-ionized gas contribute significantly for low 
strength lines. Due to our simulation resolution lines with 
EW < 20 mA could be underestimated. Following Rao & 
Turnshek (2000) and Tripp & Savage (2000) we estimate 
that at EW > 40 mA where collisional ionization domi- 
nates, the gas probed by O VI absorption lines is about 
10% of the total baryonic matter. This means that the 
O VI absorption lines can actually probe about ^ 30% of 
the WHIM gas, since the fraction of the baryons predicted 
to be in the WHIM gas is about one third of the total 
baryonic matter (Dave et al. 2001). The remaining ~ 70% 
WHIM gas is likely to be too hot to be detected via O VI 
absorption: Dave ct al. (2001) predicted that at z = the 
baryon fraction peaks at around 4 x 10^ K, in which O VII 
or O VIII should be dominated. 
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Fig. 3. — Cumulative number of O VI absorption lines per unit 
redshift vs. line equivalent width. Dashed line: model A; solid line: 
model B. Lines produced by coUis ionally ionzed gas dominate at 
EW > 40 mA. Observational data are plotted with la error bars 
(Tripp & Savage 2000; Tripp, Savage & Jenkins 2000). 

Our result is consistent with what Cen et al. (2001) 

find in their simulations, especially in that they find at 
EW > 35 mA lines from coUisionally ionized gas outnum- 
ber lines from photo-ionized gas. However, overall we pre- 
dict more O VI absorption lines. 

We find the Doppler b parameter of the O VI absorp- 
tion line can provide important clues in distinguishing 
photo- and coUisional-ionization. By fitting the spec- 
tra, the Doppler b parameters arc obtained for all the 
O VI absorption lines. We find that the distribution of 
the line widths is well approximated by a log-normal: 
/ oc exp [-log^(6/6/„)/2(Tf„]. In figure 4, we plot 6;„ 
against equivalent width cutoff for our two models. In 
the collisional ionization case (model A), the distribution 
is nearly constant with 6(„ ~ 20km/s. On the other hand, 
when we include a photo- ionizing background (model B), 
the distribution depends on the strength of the line. For 
lines with EW > 70mA, the value of bin is still around 
20 km/s, but drops quickly for smaller lines. This re- 
flects the fact that more photo-ionized lines contribute for 



^A recent observation of the low redshift quasar E 1821+643 (Tripp et al. 2001) reveals that in additional to a narrow H I Lya line, a 
relatively broaden component appears with & » 85 km s~^ in the z = 0.1212 O VI absorption systems. 
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Fig. 4. — The peak of the Doppler b distribution (fe;„) for Hnes 
with a EW larg er than the corresponding value. The value of bi„ 
is obtained by fitting each di stribution with a lognormal functional 
form (solid line: model A; dashed line: model B). The O VI ab- 
sorption lines that arc produced by collisionally ionized gas ha ve a 
mean b value of ~ 20 km a^^. After including photo-ionized gas , 
the characteristic b value drops as low as ~ 5 km s^^. 

A statistically complete sample docs not yet exist to 
compare against these predictions. Recently Tripp (2001) 
compiled the first batch of ~ 15 O VI absorption lines de- 
tected in a total of five quasars. They found that the ma- 
jority of the O VI lines have 6-valucs of ^ 22 km s~^, im- 
plying a coUisional-ionization origin. If this result contin- 
ues to hold for low equivalent width lines, it may be telling 
us something either about the strength of the photoioniz- 
ing background or the metallicity of the low-temperature 
filaments. Either case (lower background or lower metal- 
licity in the cold gas), would produce fewer photo- ionized 
lines relative to collisionally ionized. 

4. CONCLUSIONS 

In this Letter wc hixvc studied th(^ properties of O VI 
absorption lines via hydrodynamic simulation. Our main 



conclusions are as follows. 

1. The projected O VI number density distribution 
shows that the majority of O VI absorption comes 
from filaments that connect virialized intersections 
or groups of galaxies. These regions typically have 
an overdensity of ~ 5 — 100. 

2. Spectral fits show that collisional ionization 
dominates in high density, high temperature 

regions, while photo-ionization dominates in low 
density, low temperature regions. 

3. The predicted cumulative distribution of O VI 
absorption lines fits the observations remarkably 
well, given the uncertainties due to limited reso- 
lution and the assumed metallicity distribution. 
At EW > 40 mA O VI absorption lines produced 
by collisionally ionized gas start to outnumber 
lines produced by photo-ionized gas, while at 
EW > 80 mA essentially all the O VI absorption 
lines are produced by collisionally ionized gas. 



4. 



Wc find that the O VI absorption lines produced 
by collisionally ionized gas have a characteristic 



Doppler b parameter of ~ 20 — 23 km 



and 



lines due to the photo-ionized gas typically have 
narrower 6-parameters. 
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